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A Stochastic Geometry Based Approach to Modeling 
Interference Correlation in Cooperative Relay Networks 

Young Jin Chun, Simon L. Cotton, Mazen O. Hasna, and Ali Ghrayeb 


Abstract- Future wireless networks are expected to be 
a convergence of many diverse network technologies and 
architectures, such as cellular networks, wireless local area 
networks, sensor networks, and device to device commu¬ 
nications. Through cooperation between dissimilar wireless 
devices, this new combined network topology promises to 
unlock ever larger data rates and provide truly ubiquitous 
coverage for end users, as well as enabling higher spectral 
efficiency. However, it also increases the risk of co-channel 
interference and introduces the possibility of correlation in the 
aggregated interference that not only impacts the communica¬ 
tion performance, but also makes the associated mathematical 
analysis much more complex. To address this problem and 
evaluate the communication performance of cooperative relay 
networks, we adopt a stochastic geometry based approach by 
assuming that the interfering nodes are randomly distributed 
according to a Poisson point process (PPP). We also use 
a random medium access protocol to counteract the effects 
of interference correlation. Using this approach, we derive 
novel closed-form expressions for the successful transmission 
probability and local delay of a relay network with correlated 
interference. As well as this, we find the optimal transmission 
probability p that jointly maximizes the successful trans¬ 
mission probability and minimizes the local delay. Finally 
numerical results are provided to confirm that the proposed 
joint optimization strategy achieves a significant performance 
gain compared to a conventional scheme. 

Index Terms —Relay network, cooperative communication, 
stochastic geometry. Interference correlation, mean local delay. 


I. Introduction 

Cooperative relaying is an effective technique for improving 
the reliability and throughput of the traditional point-to-point 
communication. The approach was first proposed by Cover 
and El Gamal in HI and revisited in El, 0- In 0, the 
authors proposed several relaying protocols, such as amplify- 
and-forward (AF), decode-and-forward (DF), and selection 
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relaying, and evaluated the performance of these relaying pro¬ 
tocols. In 0, the authors combined relaying with space-time 
coding. More advanced relaying protocols, such as compress- 
and-forward and compute-and-forward, were introduced in 
0-IIIO1. In compress-and-forward relaying, the relay observe 
a vector-quantized signal and forward this information to 
the destination. In compute-and-forward (CF) relaying, the 
relays decode and forward linear equations of transmitted 
messages using the noisy linear combinations provided by 
the channel. Once the destination node receives enough linear 
combinations, it can successfully detect the desired messages 
m. Most of the previous work on cooperative relaying has 
focused on orthogonal channel allocation or noise-limited 
fading environments, ignoring co-channel interference. How¬ 
ever, this assumption is not realistic due to the high spectral 
reuse in practical wireless networks. The effect of co-channel 
interference on cooperative relaying has been studied in HD- 
ifTSll assuming fixed interfering node deployments. In ifTTll . the 
authors considered a relay network with interference affecting 
only the relay. While AF and DF relaying in an interference 
network were investigated in ifT^ and ifTsl . respectively. 

In a real wireless network, however, it is more practical 
to assume that the interference and node locations are random 
due to mobility. Moreover, analyzing a specific instance of the 
network with a fixed node deployment does not provide a gen¬ 
eral result. Instead, a statistical statement about the ensembles 
of all possible node deployments is much more beneficial to 
assist in understanding the network performance. In this con¬ 
text, stochastic geometry has recently gained much attention. 
This method models interference in the network by treating 
the locations of the interferer as points distributed according 
to a spatial point process m. Such an approach captures 
the topological randomness in the network geometry, provides 
well-established mathematical tools, allows high analytical 
flexibility and achieves an accurate performance evaluation. A 
common assumption in most of the related works is that the 
interfering nodes are distributed according to a homogeneous 
Poisson point process (PPP). Two important properties of PPP 
are stationarity, i.e., its distribution is invariant under arbitrary 
translation, and the Slivnyaks theorem, which means that 
conditioning on a certain point does not change the distribution 
of the process m. Due to these two properties, the PPP 
model is analytically tractable and flexible. The probability 
generating functional (PGFL) of PPP is derived in closed- 
form, and the distribution of the inter-node distance is known 
ifTSll . The Laplace transform of the interference in a PPP 
network as well as the probability density function of the 
aggregated interference were analyzed for Rayleigh fading 
channels in na, El. The outage probability and average 
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achievable rate of heterogeneous cellular networks (HetNets) 
were evaluated for PPP and Poisson cluster process (PCP) 
in llSl and ll^ . respectively. Cooperative relaying with PPP 
distributed interfering nodes has been investigated in m- 
lEi]. In na, the authors derived the throughput scaling law 
for opportunistic relay selection, while in lfT9l . decentralized 
relay selection schemes based on the location information or 
the received signal strength were proposed. To ensure a certain 
quality of service (QoS) at the destination node, the authors 
of II 20 II have defined a QoS region for random relay selection. 

Although assuming a homogeneous PPP offers a convenient 
method to model a network with uniformly distributed inter¬ 
fering nodes, using the PPP alone is not enough to accurately 
capture the real aspects of practical wireless networks. An 
important example of this occurs in the statistical behavior 
of the aggregated interference when the non-desired signals 
are correlated in either the space or time domain. Interfer¬ 
ence usually originates from a set of transmitters sharing 
common randomness, causing correlation in the aggregated 
interference. The impact of interference correlation was not 
properly reflected in the stochastic geometry based modeling 
process until very recently ll22l - ll^ . It was first addressed 
in II 22 II for random wireless networks with PPP distributed 
nodes, where the spatio-temporal correlation coefficient was 
introduced. The temporal correlation coefficient was evaluated 
in for general network models, including the static and 
random node locations for various traffic types. The diversity 
loss of a multi-antenna receiver due to interference correlation 
was analyzed in ll24l . For cooperative relaying, the interference 
correlation occurs between different receivers that are closely 
located to each other. In ||25l, the authors assumed the PPP 
model for the interfering nodes and proved that the temporal 
and spatial correlation of interference significantly degrades 
the performance of the cooperative relay. 

One effective method to reduce the interference correlation 
is to intentionally induce man-made randomness by using 
random medium access, i.e., increasing randomness in the 
MAC domain, more specifically using frequency-hopping mul¬ 
tiple access (FHMA) and ALOHA, which helps to reduce 
the effect of interference correlation. In ALOHA, each node 
transmits with a certain probability p. Decreasing the transmit 
probability increases the uncertainty in the active interfering 
nodes and reduces the interference, thereby reducing corre¬ 
lation. In 1261, the authors analyzed the local delay, which 
is the time it takes for a node to successfully transmit to 
a nearby neighbor, using FHMA and ALOHA on PPP and 
determined the optimal number of sub-bands in FHMA and 
the optimal transmit probability in ALOHA that minimizes the 
local delay. Stamatiou and Haenggi ll27l applied this approach 
to a multi-hop relay network and evaluated the local delay of 
time division multiple access (TDMA) and ALOHA protocols. 
However, the authors assumed the nodes to be aligned on a 
one dimensional straight line which limits the application of 
their work. 

To the best of our knowledge, with the exception of lIZTl . 
there is no previous work that has considered the effect of 
the MAC protocol on interference correlation for a cooper¬ 
ative relay network. This motivates us to consider a more 


realistic relaying network where the nodes are spread over 
a two dimensional vector space and to evaluate the effect of 
interference correlation on the communication performance. 
In particular, for the first time, we compare the successful 
transmission probability and local delay of a relay network, 
which are subject to both correlated and uncorrelated interfer¬ 
ence. We also propose an efficient optimization strategy that 
jointly maximizes the success probability and minimizes the 
local delay. Specifically, we provide the following theoretical 
contributions. 

1) We analyze the successful transmission probability and 
local delay of a relay network with PPP interfering 
nodes. We consider both interference correlated and 
uncorrelated cases. 

2) We determine the necessary and sufficient condition to 
jointly maximize the success probability and minimize 
the local delay, and derive the optimal transmission 
probability p that achieves this optimality. 

3) We propose an optimization strategy that numerically 
finds the optimal p, and then compare the computational 
complexity of our proposed strategy to that of the 
conventional brute-force method. 

4) We provide numerical results to validate the analysis, 
evaluate the performance of the interference correlated 
and uncorrelated cases, and compare the performance 
gain that our proposed optimization strategy achieves 
compared to the conventional scheme. 

The rest of this paper is organized as follows. In Section II, 
we describe the system model and introduce the mathemati¬ 
cal background for the interference analysis performed here. 
We derive the successful transmission probability and local 
delay for the interference correlated and uncorrelated cases in 
Section III. In Section IV, we determine the necessary and 
sufficient condition to achieve the joint optimality between 
the success probability and local delay, and propose an opti¬ 
mization strategy based on iterative numerical search. Section 
V provides some numerical results based on our approach. 
Finally, Section VI concludes the paper with some closing 
remarks. 

H. System AND Mathematical Models 
A. System Model 

We consider a three node relay network, consisting of a 
source, relay, and destination, where multiple interfering nodes 
simultaneously transmit during each time slot, as illustrated 
in Fig. [T] We denote the source, relay, destination, and the 
interfering node as s, r, d, and x, respectively, where the 
notations denote both the nodes and their coordinates. We 
assume that the interfering nodes are randomly distributed 
according to a Poisson point process (PPP) $ of intensity 
A, i.e., X G $, and the destination is located at the origin 
d = (0,0). The source, relay, and interfering nodes transmit 
with power Pg, Pr, and P^, respectively. 

We use the ALOHA protocol with transmit probability p on 
each time slot, i.e., the source attempts to access a slot with 
success probability p. Given the channel access, the source 
transmits its packet to the relay. We follow the approach of 
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Fig. 1: System Model 


II 25 I and by assuming that a certain time is reserved for 
the relay transmission right after the source transmission, so 
that the relay forwards the packet to the destination within the 
same time slot. If the relayed transmission fails, the source 
attempts to re-transmit during the next time slot. We assume 
that the link between source and destination is unreliable, so 
the transmission occurs only through the relay, i.e., there is no 
direct link between s —^ d. 

The distance between arbitrary node i and j is denoted 
by II* ~ ill and the path loss function between two nodes 
is given by l{i,j) = ||i — j||““, where a > 2 is the path 
loss exponent. The channel links are assumed to be subject to 
independent and identically distributed (i.i.d.) Rayleigh fading 
with mean one, where the channel coefficient between node i 
and j is denoted by hij. The additive noise w is assumed to 
be complex Gaussian distributed with mean zero and power 
spectral density Nq. 


B. Mathematical Model 

Let denote the set of active interfering nodes in time slot 
k, i.e., € $. The aggregated interference at the destination 

during time slot k is 

^k,d — Px ^ ^ kxd^ix, d?j\{x G 

xg$\{s} 

and the interference at the relay during time k is given by 


^k,r — Px ^ ^ kxr^{x,T^\{x € 


( 2 ) 


where 1( ) is the indicator function. Then, the signal to 
interference plus noise ratio (SINR) between s ^ r and r ^ d 
links during time slot k is 


^k,sr 


Pshsrl{s,r) 

No + Ik,r 


^k.rd 


d) 

Nq + Ik,d 


Pshsrl{s,r) 

N + Il^ 


Prkrd^(,’l^, d) 

N + Iid ’ 


(3) 


where 4 p^/p^, p^ A ^ A j' A jjp^^ 

In this paper, we consider the cases of both correlated 
and uncorrelated interference. For correlated interference, we 
assume that the interference at the destination and the relay 
originate from the same set of interfering nodes, i.e., <i>fc in 
O and (l2]i. For uncorrelated interference, the interference at 
the destination and the relay are generated by two different 
sets of interfering nodes, i.e., <I>fe ^ such that the inter¬ 
ference model in ([Til, dU must be modified. The uncorrelated 
interference case is analyzed in detail in Section III-C. 

III. Analysis of the System Measures 
A. Packet Delivery Probability 

The packet delivery probability, i.e., the successful trans¬ 
mission probability between s ^ r ^ d, is now derived in 
Theorem 1. 

Theorem 1. The packet delivery probability P{C) of the three 
node relay network which is subject to correlated interference 
is given by 


P{C) = pexp{—Xp'tp{r) — B), (4) 

where ip{r) is defined in (|7|), B = N “ 

9/l{i,j) = e\\i-j\\°‘. 

Proof: See Appendix I. ■ 


B. Mean Local Delay 

The source attempts to re-transmit if the destination fails to 
receive the packet during the previous time slot. In general, a 
successful transmission during each time slot is a dependent 
event due to the correlated interference ll2^ . However, if we 
consider the conditional success event for a given <I>, the 
randomness stems only from the channel fading coefficient 
and the ALOHA protocol which are independent variables for 
each time slot. Therefore, the success event in different time 
slots given $ are independent with probability P (C$) in (l33]) . 

Let us define the local delay as the number of time slots 
required until a successful transmission occurs. Then, the 
local delay given $, denoted by A$, can be represented as 
a geometric random variable written below 

P = k) = - P P {C^), k>l. (5) 

The mean local delay D{p) = E[E(A$)] averaged over all 
possible $ is now derived in Theorem 2 below. 

Theorem 2. The mean local delay D{p) of the three node 
relay network for correlated interference is given by 

D{p) = ^ exp {Xpip{r) + B ), (6) 

where (p(r) is defined in (|S|) and B = N 

Proof: See Appendix II. ■ 
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V^(r) = f 

JR^ 


1 - 


Ssrljx.r) ^ _|_ er.dl{x,d) '^ 


dx = 


I R'i 


1 - 


-(- g«>-|k-r|h° ^ ^2^ _|_ grdikl 


V{r) ^ [ - 
Jr'^ ^ 


fix) 


^2 1 + (1 -p)/(x) 


dx, f{x) = 1 + 


A / -, , Ssrlix,r)\ f OrdKXyd) 


1 + 


Pr 


- 1 , 


dx, (7) 


( 8 ) 


f^uir) = [ 
Jr? 


1 - 


/i I 0srl{x,r) 
V Ps 


dx = 


1 


‘Puir) = [ ( \ dx, g{x) = 

Jr 2 1 + {l-p)g{x) 


R? j^||x — r||“ + 1 
ds7-/(x, T ^ 


-dx, 


P, 


(9) 


( 10 ) 


Ep 

g“ (^p-ffc,r + -^ A,d) 

= n ® 

esrlG,r) 

6 

n ® 




a:e<I>\{s} 


a:e$\{s} 

L J 


= n 

xG$\{s} 


P 

1 I 

Ps 


+ 1- p 


n 

xe$\{s} 


P 

1 _j_ ^rddX' idj 

-Pr 


+ 1-P , 


( 11 ) 


C. Results for Uncorrelated Interference 
For uncorrelated interference, the interference at the des¬ 
tination and the relay are generated by two different sets of 
nodes: the aggregated interference at d and r are given by 

Jk,d — Px ^ ^ dxd^ixt d)\{x G ^fc), 

x6$\{4 ^ ^J2) 

Ik,r = Px ^ hxrl{x,r)l{x G 

x£i\{s} 

where $fc and are two different sets of active interfering 
nodes, i.e., Then, the expectation in the last equality 

of (133b is evaluated in (HUl and the corresponding packet 
delivery probability and the mean local delay for independent 
interference are derived in Theorem 3 below. 


Theorem 3. The packet delivery probability P{C) and the 
mean local delay D{p) of the three node relay network for 
uncorrelated interference are given by 

P((7) =pexp - A7rpC’(d) ~ 

D{p) = i exp (^Ap(p„(r) -f , (14) 

where '0u(r), Puix), and B are defined in ( 170b . and 
Theorem 1, respectively. 

Proof: See Appendix III. ■ 


Remark 1. Without a MAC coordination, i.e., p = \, the mean 
local delay becomes infinity as indicated in 0261 . We observe 
the same result on cooperative relay networks. By using (|6]l, 
the mean local delay for correlated interference with p = 1 
case is given by 


Dip) = exp (a p(r)|p^i -f p) . 


Pix)\p=i is lower bounded as follows 


Pix)\p=i = [ fix)dx > ^ / 


IR^ 
drd 
Pr 


P, 

-“dx = 


/(x, d)dx 
‘27T9rd 


Ir'^ 


Pr 


oo (16) 

r^-“dr. 


' 7 


' r—0 


where we applied (0 in the first equality and transformed the 
Cartesian coordinates to Polar coordinates, i.e., x —>■ re™, in 
the last equality. Then, ip{r) diverges to infinity at p = 1. 
Similarly, for uncorrelated interference, the term within the 
exponential of (fT4l i diverges to infinity. Hence, the mean local 
delay of the three node relay network becomes infinity at p = 
1 . 


TV. Optimization Strategy 

In this section, we determine the ALOHA transmission 
probability p that optimizes the system measure in Section 
III. We determine the necessary and sufficient condition for 
optimality, propose an iterative method to find the optimal p* 
that meets both conditions, and compare the computational 
complexity of our proposed optimization strategy to that of 
the brute-force search method. 


A. Minimizing the Delay 

We minimize the mean local delay 79 (p) as follows 

Find p* that minimizes D(p), i.e., —= 0. (17) 

o<p<i up 

Corollary 1. The optimal transmission probability p* that 
minimizes 79 (p) is the solution of the following condition 

1 dtp{r) 

- = r{r) + v-g^ 


(15) 


(18) 
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for correlated interference, and 


for uncorrelated interference. 


'kC(6)(1-pS) 



Proof: For brevity, we only prove the interference cor¬ 
related case, however the proof for uncorrelated interference 
can be found using a similar approach. By using Theorem 2, 
the first derivative of D{p) is obtained as 


dD{p) 

dp 


D{P) 

P 


1 — Xp 


(^ip{r) +p 


dp{r) \ 

dp J 


( 20 ) 


Since D(p) has a positive value, the optimal p* achieves (fT^ . 
The second derivative of D{p) at p* is given by 


d^Djp) _ D{pl 

Qp2 p2 

p* 


1 + Xp^ 



+ p 


af^£fr)\ 
dp^ J 


( 21 ) 


where we applied (fT^ to simplify the expression. The deriva¬ 
tives of (p{r) are given by 


9p(r) ^ r fjx)^ 

dp 7^. (l + (l-p)/(x))2 ’ 

d'^pjr) ^ r f(^x)^ 

dp^ 7^, (l + (l_p)/(x))3 

Since f{x) is a non-negative function on x, and ^ 

are both non-negative. Hence, ^ in (|2TI) has a positive 
value at p* and the transmission probability p* minimizes the 
mean local delay D{p). This completes the proof. ■ 


B. Joint Optimization of the Throughput and Delay 

Let us define the utility function U{p) = as the 

ratio of the network throughput pP{C) to the mean local 
delay D{p). By maximizing U, we can jointly maximize the 
throughput and minimize the delay as follows 

dU 

Find p* that maximizes { 7 , i.e., = 0. (23) 

o<p<i dp 

Corollary 2. The optimal probability p* that maximizes the 
utility U is the solution of the following condition 


— = ip{r) +p{r) +p 
Xp 

for correlated interference, and 


d(p(r) 

dp 


3 

Xp 


= tpu{r) + Pu{r) +p 


dpujr) 

dp 


+ nC{5) 




{l-p5) 

(l-p)2-^ 


(24) 


(25) 


for uncorrelated interference. 


Proof: Again, we only prove the interference correlated 
case, since the uncorrelated interference case can be obtained 
in a similar manner. By using Theorem 1 and 2, the first 
derivative of U can be expressed as follows 


dU(p) _ U(p) 

dp p 


3 - Ap ( tf{r) + ip{r) + p^^ 


(26) 


Since U{p) has a positive value for p > 0, the optimal p* 
achieves (l24l i. The second derivative of U{p) at p* is given by 


dp'^ 


u(p) 

P'^ 


3 -|- Xp^ 


^dpfr) 


P dp^ 


(27) 


where we applied (l24li to simplify the expression. Since 

and are both non-negative, ^ in (l27l i has a negative 

value at p* and the transmission probability p* maximizes the 
utility U{p). This completes the proof. ■ 


C. Proposed Optimization Strategy 


Algorithm 1 Find p* that achieves (fTsTl. (fT^. (l24ll. 


Require: 0 < po < 1, eo > 0 


procedure Newton-Raphson Iteration 
P<^Po 

while |A(p)| > eo do 

A(p) 

P^P-A^ 

end while 
return p* p 
end procedure 


The optimal conditions in Corollary 1 and 2 are implicit 
functions of p: the effect of the interference and the ALOHA 
transmission probability p are coupled inside the integral term 
of (p(r) and Pii(r), which makes it hard to explicitly calculate 
p* from (fTSl) . ( fT^ . (l24l i. and (l25l) . Hence, we adopt an iterative 
search method to numerically find the optimal p*. Specifically, 
we use Newton-Raphson iteration as outlined in Algorithm 1 
1^ . First, we denote A(p) for (l24l i as 

A(p) = ^-'!/’(r) - p(r)(28) 

Similarly, A(p) for (fTSl l. ( fT^ , and (l25l) are defined by subtract¬ 
ing the terms on the left-hand side to that on the right-hand 
side. The derivative of A(p) for (l24l i is given by 


a'(p) 


3 ^dip{r) d‘^ip{r) 
Ap^ dp ^ dp"^ 


(29) 


where the derivatives of p(r) are derived in (l22l i. 

The iterative search begins by setting the initial transmission 
probability 0 < po < 1 and the termination criteria cq > 0. 
The ALOHA transmission probability is updated as 


Pm-\-l — Pm 


■^{Pm} 
^'{Pm) ’ 


m : iteration index, 


(30) 


until A{pm) becomes arbitrarily close to zero. If the trans¬ 
mission probability pm satisfies the condition |A(pm )l < Co, 
the iterative search stops and returns p* = p^ as the optimal 
transmission probability. 

Conventional approaches 1291 for numerically finding p* 
include the brute-force method where the utility C/(p) is 
evaluated for all possible 0 < p < 1 and compared to 
find the optimal value and the quick-sort method where each 
adjacent numbers U{pi) and U{p 2 ) are compared and sorted 
until the maximum value is reached. The ratio between the 


p 
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computational complexity of the Newton-Raphson iteration 
and that of the brute-force method is O{\ogn)/0{n^) for n 
digit precision, whereas the ratio between the computational 
complexity of the Newton-Raphson iteration and that of the 
quick-sort method is O(log n)/0(n log n) 1^ . Hence, the op¬ 
timization strategy proposed here is computationally efficient 
and effective compared to the conventional schemes, such as 
brute-force method or quick-sort method. 


V. Numerical Results 


In this section, we present numerical examples through 
which we compare the performance of the interference cor¬ 
related and uncorrelated cases. Fig. |2] shows the success 
probability, mean local delay, and utility calculated using (|4|i, 
©, O, (O, and U(p) = and evaluated over different 

relay locations. We use a three node relay network where the 
nodes are located at s = (2,0), d = (0,0), r = {rx,0). The 
interfering nodes are randomly deployed by using PPP with 
node density A. We assume that the nodes transmit with power 
Ps = 5 dB, Pr = 5 dB, Px = 5 dB, the noise power is Nq = 1, 
the path loss exponent is a = 4, and the SINR threshold is 
0=1. The dotted and solid curves are the numerical results for 
the interference correlated (IC) and interference uncorrelated 
(lU) cases, respectively. It is clear that the IC achieves a better 
performance than the lU in terms of success probability, mean 
local delay, and utility. For a small interfering node density 
A, the performance gap between the IC and lU cases is very 
small, however, as the node density increases, the gap becomes 
more evident. We also note that for a small node density, using 
a high transmission probability achieves a better performance 
in terms of the success probability, mean local delay, and 
utility. For a large node density, more nodes will interfere with 
each other, increasing the mean local delay and decreasing the 
success probability. So, it is optimal to use a low transmission 
probability for large node density case. 

Fig. [3 shows the success probability versus node density, 
transmission probability, and SINR threshold where we fixed 
the relay location at the middle, i.e., r = (1,0). The dotted and 
solid curves are again the numerical results for the IC and lU 
cases with a fixed transmission probability p. For the curves 
without a line and only the markers, we determined the optimal 
p* that maximizes the utility U {p) using Corollary 2 for each 
node density A and SINR threshold 9. Then, we evaluated the 
corresponding success probability for the given (p*. A, 9). It is 
clear that the IC case achieves better performance than the lU 
case and the performance gap increases as the node density 
increases. By using similar method that was used in Theorem 
1 and 3, the link success probability can be derived as m 


p [lk,Td > 9] = exp 


P i'jk.sr > 9]= exp 




(31) 


which are decreasing function of the SINR threshold 9. So, 
each link falls in outage with a high probability for large 9 
and the gap between IC and lU becomes increasingly wider 


as the node density increases. Since we jointly optimized the 
network throughput pP(C) and the mean local delay D{p), the 
optimized p do not achieve the maximum success probability 
for all cases. However, our proposed optimization strategy 
still achieves a significant performance gain, compared to the 
conventional fixed transmission probability case. 

Similarly, Fig. 0] shows the mean local delay versus node 
density, transmission probability, and SINR threshold where 
the relay location was fixed at the middle, i.e., r = (1,0). The 
IC case achieves a lower mean local delay than the lU case 
in most of the scenarios and the performance gap increases 
for a large node density. Fig. 4(b) illustrates the convexity 
of the mean local delay, i.e., as the node density increases, 
the transmission probability that minimizes the mean local 
delay decreases. Again, although the optimized p determined 
by Corollary 2 does not achieve the minimum mean local delay 
for all cases, it still achieves a significant performance gain 
compared to the conventional fixed transmission probability 
case. 


VI. Conclusion 

In this paper, we have considered a three node relay 
network subject to interference generated by nodes which were 
distributed according to a PPP. Furthermore, we evaluated 
the packet delivery probability and the mean local delay for 
both interference correlated and uncorrelated cases. Based on 
these analytical derivations, we determined the necessary and 
sufficient condition to jointly maximize the packet delivery 
probability and minimize the mean local delay. The trans¬ 
mission probability p that achieves this optimality was also 
found. Specifically, we observed that correlated interference 
achieves better performance than uncorrelated interference in 
terms of success probability, mean local delay, and utility. For 
a small node density, the performance gap between IC and lU 
is very small, however, as the node density increases, the gap 
becomes more evident. Also, for a small node density, using 
a high transmission probability p achieves better performance. 
However for a large node density the opposite is true, i.e., a 
low transmission probability should be used. Based on these 
observations, we have proposed a computationally efficient 
optimization strategy that finds the optimal p using an iterative 
search. Finally, we have also provided numerical results to 
prove the performance gain that our proposed optimization 
strategy achieves compared to a conventional scheme. 


Appendix I 

In this appendix, we prove Theorem 1. Let denote the 
successful transmission event between s —>■ r —>■ d conditioned 
upon the PPP $. Using Q and a pre-defined SINR threshold 
9, C <5 can be written as 

— {^k,sr ^ 9 ^k,rd ^ 


(32) 
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and the conditional success probability given $ is 

P(C$) = pP ['yk,sr > On ')k,rd > 

= pP Pshsr > dsr{N + 

n PrKd > 0,d(iv + /'^1$ 

= pEj^ 

iCr ( ^sr I_ 

Ps + Pr JE 




(33) 


= pe 


_ ^sr J _ ^rd 

g Ps Pr 


Ej/ 


= n ® 

a:e$\{s} 

= n 


— {t' I ^rd / 

- + ftxd) l(a:e$fc) 


_j_ ^sr^(^ 5 ^) ^ I ^rd^i.^t^') 


+ 1 -P 


P(C) =E[P(C$)] 

_jv( 

= pe ( 'P'i J’r / E 


where 


j{x) = 




n 

xG$\{s} 


+ 1- p. 


G[r;(a:)] = E 


n 


.a:€$ 


= exp J [3 ~ u(x)](ix^ , 


Appendix II 

In this appendix, we prove Theorem 2. The mean of the 
geometric random variable in (|5]l is E [A$] = l/P(C<i>). Then, 
the mean local delay for all available $ is given by 

D(p) 4E|1(A,)| =E ' 


P(C.)J ■ 

By substituting ([33]) and ( l34l i into the last equality in ( [38] ). 
applying the PGFL of PPP ( l37b . the mean local delay can be 
expressed as follows 


where the first equality is due to the ALOHA transmission 
probability p and the distribution of hij, i.e., P{hij > x) = 
exp(—x), is applied to the third equality. 

For correlated interferenc^H the interference is generated 
from the same set By substituting dD and (|2]) in the last 
equality of (i33T l and averaging it with respect to the fading 
coefficient hxi and the ALOHA protocol the expectation 
in the last term of (i33T l can be expressed as 


Dip) = —E 
P 

1 


n 


1 


vix) 


(39) 


= - exp P — A j 
P \ Jr^ I 


1 - 


1 


v(x) 


dx 


Since the integral in ( l39b can be simplified as 

' dx, 


-pp(r) = 1 - 

Jr^ L 


’{x)_ 


(40) 


(34) 


we obtain the mean local delay in (|6]) by substituting (l40l i into 
(i39] l. This completes the proof. 

Appendix HI 

In this appendix, we prove Theorem 3. The conditional 
success probability given $ is 


Using (i33T l and (l34l i. the packet delivery probability, i.e., the 
average success probability over all possible $, can now be 
expressed as follows 


P(C< 5 ) = pe ^Ej' 

n 


(41) 


= pe 


dl(.x) P2ix), 

where the first equality follows by (i33T l, the second equality 
follows by (im i. pi(x) and pL 2 ix) are denoted by 

P 


piix) = 


(35) 


1 + 


0arl{x,r) 


1 -P, 


P2ix) = 


P 


1 + 


0rdl{x,d) 


+ 1 - P. 


(42) 


(36) 


The expectation term in (l35l l is referred to as the probability 
generating functional (PGFL) and can be evaluated for PPP as 

(na 


By averaging (EB over all possible $, the packet delivery 
probability for uncorrelated interference can be expressed as 


(43) 


PiC)=pe-^E 

n 

E 

n /^2(a;) 


2 :e<[>\{s} 




= pexp ~ ^ 

X exp ( —A 


[1 — idi{x)]dx 


(37) 


where the integral in (ITTT i is a two dimensional integral since 
X is a coordinate. Hence, we obtain the packet delivery 
probability in (|4]) by substituting (13^ and (l37l i into dTSl l and 
using the notation ipir) in (j?]). This completes the proof. 

* In Section III-A and B, we only consider correlated interference. The 
results for uncon'elated inference are summaiized in Section III-C. 


7i?2 

/ [l-/i 2 (x)]dx 

V 

where the first equality follows by (1351) and the PGFL of PPP 
( [37] ) is used for the second equality. The first integral term 
in (EB can be represented as /^ 2[1 — ^i(x)]dx = ptpuir) by 
using straight forward calculus and the second integral term 
can be expressed in closed form as follows 

1 


/ [1 - fi2ix)]dx = p — 

JR^ 1 


-dx 


0rd 


= 27rp 


p=0 1 + 


-dp = TvpCiS) ( 


( 44 ) 


















































where we applied the notation of ^ 2 {x) in (1421) to the first 
equality, transformed the Cartesian coordinates to Polar coor¬ 
dinates, i.e., X pe“, in the second equality, and applied the 
following integration equality in the last inequality ll^ 




1 -f qx’" 


-dx = 




q-^c(f^), C(5) = ^i—. (45) 

\iy/ sinc(o) 


(fUt follows by substituting (l44l i and (|9]l into ( l43l l. 

Similarly, the mean local delay can be derived as follows 




D{p) = —E 
P 


n 


Plix) 


E 


= - exp i? — A 
P 


X exp I —A 




Jr^ 
1 - 


1 - 


n 

xG'I’Xfs} 

dx 


P2ix) 


P2{x) 


Hi{x) 

dx 


(46) 


where the two integral terms in ( l46l l can be evaluated as 


f 

'_R2 

1 

ri{x)_ 

f 

1 

P2{x)_ 


dx = -p(pu{r), 


^ _ TTpC(S) 
(1-py-s 



(47) 


by using the same integration technique in (l44l) and (|45] |. 
Hence, we obtain the mean local delay in (fT4l) by substituting 
Wh and ([Tol l into (l46l l. This completes the proof. 
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(a) (b) (c) 

Fig. 2: (a) Success probability, (b) mean local delay, and (c) Utility versus relay location. 



(a) (b) (c) 


Fig. 3: Success probability versus (a) node density, (b) transmission probability, and (c) SINR threshold. 



Fig. 4: Mean local delay versus (a) node density, (b) transmission probability, and (c) SINR threshold. 










































































































